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Fluctuation Cumulant Behavior for the Field-Pulse Induced 
Magnetisation-Reversal Transition in Ising Models 
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The universality class of the dynamic magnetisation-reversal transition, induced by a competing 
field pulse, in an Ising model on a square lattice, below its static ordering temperature, is studied 
here using Monte Carlo simulations. Fourth order cumulant of the order parameter distribution 
is studied for different system sizes around the phase boundary region. The crossing point of the 
cumulant (for different system sizes) gives the transition point and the value of the cumulant at the 
transition point indicates the universality class of the transition. The cumulant value at the crossing 
point for low temperature and pulse width range is observed to be significantly less than that for 
the static transition in the same two-dimensional Ising model. The finite size scaling behaviour in 
this range also indicates a higher correlation length exponent value. For higher temperature and 
pulse width range, the transition seems to fall in a mean-field like universality class. 

PACS numbers: 05.50.+q; 05.70.Fh 
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I. INTRODUCTION 

The response of a pure magnetic system to time- 
dependent external magnetic fields has been of current 
interest in statistical physics 0, 0, Q These studies, 
having close applications in recording and switching in- 
dustry, have also got considerable practical importance. 
These spin systems, driven by time-dependent external 
magnetic fields, have basically got a competition between 
two time scales: the time-period of the driving field and 
the relaxation time of the driven system. This gives rise 
to interesting non-equilibrium phenomena. Tome and 
Oliveira first made a mean-field study jj| of kinetic Ising 
systems under oscillating field. The existence of the dy- 
namic phase transition for such a system and its nature 
have been thoroughly studied using extensive Monte- 
Carlo simulations. Later, investigations were extended 
to the dynamic response of (ferromagnetic) pure Ising 
systems under magnetic fields of finite-time duration [5( . 
All the studies with pulsed field were made below Tj?, the 
static critical temperature (without any field) , where the 
system gets ordered. A 'positive' pulse is one which is 
applied along the direction of prevalent order, while the 
'negative' one is applied opposite to that. The results for 
the positive pulse case did not involve any new thermo- 
dynamic scale 5J. In the negativepulse case, however, 
interesting features were observed j5|: the negative field 
pulse competes with the existing order, and the system 
makes a transition from one ordered state characterised 
by an equilibrium magnetisation +mo (say) to the other 
equivalent ordered state with equilibrium magnetisation 
—mo, depending on the temperature T, field strength 
hp and its duration At. This transion is well studied 



in the limit At — > oo for any non-zero value of h p at 
any T < T®. This transition, for the general cases of fi- 
nite At, is called here 'magnetisation-reversal' transition. 
Some aspects of this transition has been recently studied 
extensively 0,0- 



II. MODEL AND THE TRANSITION 

The model studied here is the Ising model with nearest- 
neighbour interaction under a time-dependent external 
magnetic field. This is described by the Hamiltonian: 



(1) 



where Jij is the cooperative interaction between the spins 
at site i and j respectively, and each nearest-neighbour 
pair is denoted by {...}. We consider a square lattice. 
The static critical temperature is T° = 2/ln(l + v2) — 
2.269... (in units of J/K B ). At T < T c °, an exter- 
nal field pulse is applied, after the system is brought 
to equilibrium characterised by an equilibrium magneti- 
sation mo(T). The spatially uniform field has a time- 
dependence as follows: 



h(t) 



-hp t < t < t + At 
otherwise. 



(2) 
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Typical time-dependent (response) magnetisation m(t) 
(= < Si >, where < ... > denotes the thermodynamic 
'ensemble' average) of the system under different mag- 
netic field h(t) are indicated in the Fig. 1. The time to 
at which the pulse is applied is chosen such that the sys- 
tem reaches its equilibrium at T (< T% ). As soon as the 
field is applied, the magnetisation m(t) starts decreasing, 
continues until time t + At when the field is withdrawn. 
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The system relaxes eventually to one of the two equlib- 
rium states (with magnetisation — m Q or +m Q ). At a 
particular temperature T, for appropriate combinations 
of hp and At, a magnetisation-reversal transition occurs, 
when the magnetisation of the system switches from one 
state of equilibrium magnetisation mo to the other with 
magnetisation —mo. This reversal phenomena at T < 
is simple and well studied for At — > oo for any non-zero 
hp. We study here the dynamics for finite At values. It 
appears that generally h p — > oo as At — » and h p — > 
as At — ► oo for any such dynamic magnetisation-reversal 
transition phase boundary at any temperature T (< Tj?). 
In fact, a simple application of the domain nucleation the- 
ory gives hplnAt = constant along the phase boundary, 
where the constant changes by a factor l/(d+ 1), where 
d denotes the lattice dimension, as the boundary changes 
from single to multi-domain region 
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FIG. 1. Typical time variation of the response magnetisation m(t) 
for two different field pulses h(t) with same Ai for an Ising system 
at a fixed temperature T. The magnetisation-reversal here occurs 
due to increased pulse strength, keeping their width At same. The 
transition can also be brought about by increasing At, keeping 
hp fixed. The inset indicates the typical phase boundaries (where 
the field withdrawal-time magnetisation m w = 0) for two different 
temperatures (sequential updating; note that for random updating 
the phase boundaries shift upwards). 

A mean field study of the problem gives a qualitative 
understanding of the diverging time and length scales 
developed near the transition boundary (in the h p - At 
plane at a fixed T < T c °). Mean-field approximation for 
the dynamics gives the equation of motion for the average 
magnetisation mi as 



transition point, gives, for uniform magnetisation, 



m(t) 



hp ( hp 

~ AT ~ VAT 



mo 



AT . 
"X]> i — (f - t ) 



(4) 



as a solution of eqn. iJSJ, for to < t < to + At. Here 
AT = T™f - T, where T c m/ = J(q = 0) is the static 
critical temperature in the mean-field approximation and 
J(q) is the Fourier transform of the interaction Jy . Due 
to application of the field h p , m(t) decreases in magnitude 
from m(to) = mo to m(to + At) = m w at the time of 
withdrawal of the pulse. Due to absence of fluctuation 
here, magnetisation relaxes back to its original value mo 
if m w is positive, or to a value —mo if m w is negative. In 
the t > to + At regime, where h(t) = 0, the magnetisation 
(starting from m w at t — t + At) relaxes back to its final 
equilibrium value ±mo, with a relaxation time 0-0 
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(T c m/ - T) 



In 



Too 



(5) 



It diverges at the magnetisation-reversal transition 
boundary, where m w vanishes. The prefactor gives the 
divergence of r at the static mean field transition temper- 
ature, and is responsible for critical slowing down phe- 
nomena at the static transition point (h = 0). The other 
factor gives the diverging time scale, at any temperature 
below the static transition temperature, where magneti- 
sation reversal occurs or m w vanishes due to appropriate 
combination of h^&nd At. The solution of the suscepti- 
bility x(q) gives 



X(q) ~ exp i 



(6) 



where the correlation length is given by 



(T c m/ - T) 



In 



toq 

m w 



(7) 



Here, too, the prefactor in \ gives the usual divergence 
at T™-> , while the other factor gives the divergence at 
the magnetisation-reversal transition point. Incorporat- 
ing fluctuations, extensive Monte-Carlo simulation stud- 
ies have also convincingly demonstrated that the fluc- 
tuation in the order parameter \m w \ and in the internal 
energy of the system grows with the system size and di- 
verges at the magnetisation-reversal transition boundary, 
where m„, vanishes. 



III. MONTE-CARLO STUDY AND THE 
RESULTS 



drrii 



tanh 



T 



Ht) \ 



(3) 



This equation, linearised near the magnetisation-reversal 



Here the Monte-Carlo study has been carried out in 
two-dimensions (square lattice) with periodic boundary 
conditions. Spins are updated following Glauber dynam- 
ics. The updating rule employed here are both sequen- 
tial as well as random. In sequential updating rule one 
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Monte-Carlo step consist of a complete scan of the lat- 
tice in a sequential manner; while in random updating a 
Monte-Carlo step is defined by N (= L 2 ) random up- 
dates on the lattice, where N is the total number of spins 
in a lattice of linear size L. Studies have been carried 
out at temperatures below the static critical tempera- 
ture (T° ~ 2.27). The system is allowed to evolve from 
an initial state of perfect order to its equilibrium state at 
temperature T. The time to is chosen to be much larger 
than the static relaxation time at that T, so that the 
system reaches an equilibrium state with magnetization 
+mo(T) before the external magnetic field is applied at 
time t = U). The field pulse of strength — h, p is applied for 
duration At (measured in Monte Carlo steps or MCS). 
The magnetisation starts decreasing from its equilibrium 
value too- The average value of the magnetisation m w 
at the time of withdrawal of pulse is noted. The phase 
boundary of this dynamic transition is defined by ap- 
propriate combination of h p and At that produces the 
magnetisation reversal by making m(to -I- At) = m w = 
from a value m(to) = mo, i.e, m w changes sign across the 
phase boundary. The phase boundary changes with T. 
The bahavior of different thermodynamic quantities are 
studied across the phase boundary. These quantities are 
averaged over 1000 — 20000 different initial configurations 
of the system. The fluctuations over the average value 
are also noted. 

Here we study the behavior of the reduced fourth or- 
der cumulant U near the magnetisation reversal tran- 
sition. This is defined as 



where (m* ) is the ensemble average of m^. (m^,) is 
similarly defined. The cumulant U here behaves some- 
what differently compared to that in static and other 
transitions: Deep inside the ordered phase m w ~ 1 and 
U -^2/3. For other (say, static) transitions the order pa- 
rameter (m w ) goes to zero with a Gaussian fluctuation 
above the transition point, giving U — > there. Here, 
however, due to the presence of the pulsed field, \m w \ 
is non-zero on both sides of the magnetisation-reversal 
transition. Hence U drops to zero at a point near the 
transition and grows again after it. 

The universality class of the dynamic transition in 
Ising model under oscillating field has been studied ex- 
tensively by investigating the critical point and the cu- 
mulant value U* at the critical point, where the cumulant 
curves cross for different system sizes (L). In that case, 
of course, the variation of U (at any fixed L) is similar to 
that in the static Ising transitions (U = 2/3 well inside 
the ordered phase and U — > well within the disordered 
phase). In fact, U* value in this oscillatory field case was 
found to be the same as that in the static case, indicating 
the same universality class 0. We observe different ba- 
haviour in the field pulse induced magnetisation-reversal 
transition case. 



We observe two kinds of distinct behavior of the 
cumulant U. Typically, for low temperature and low 
pulse-duration region (see the inset in Fig. 1) of the 
magnetisation-reversal phase boundary, the cumulant 
crossing for different system sizes (L) occur at U* ~ 0.42 
to 0.46 (see Fig. 2). As mentioned already, we have 
checked these results for both sequential and random up- 
dating. Specifically, for T = 0.5 and At = 5, (see Fig. 
2c) we find the transition point value of h p ~ 2.6, to 
be smaller than the value (~ 1.9) for sequential updat- 
ing. However, the value of U* at this transition point is 
again very close to about 0.44. This indicates that updat- 
ing rule does not affect the universality class (U* value), 
as long as the proper region of the phase boundary is 
considered. For relatively higher temperature and pulse- 
duration region of the phase boundary, the crossing of 
U for different L values occur for U* ~ 0+. This is 
true for both sequential (Fig. 3a, b) and random (Fig. 
3c) updating. It may be noted that the phase boundary 
changes with the updating rule, as the system relaxation 
time (which matches with the pulse width at the phase 
boundary) is different for sequential and random updat- 
ing 0. 
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FIG. 2. Behavior of U near the transition, driven by (a) T at a 
fixed value of h p (=1.9) and At (=5) with sequential updating, 
(b) hp at a fixed value of T (=0.5) and At (=5) with sequential 
updating, and (c) h p at a fixed value of T (=0.5) and At (=5) 
with random updating, for different L, averaged over 1000 to 20000 
initial configurations. The fluctuations are smaller than the symbol 
size. The insets show the typical behavior of the magnetisation m w 
at the time of withdrawal of the field pulse by varying (a) T at a 
fixed h p and At, for L = 100 and 800, (b) h p at a fixed T and At, 
for L = 100 and 400, (c) h v at a fixed T and At, for L = 50 and 
200; m w = at the effective transition point, (d) Finite size scaling 
study in this parameter range: the effective T c or hp 1 values (see the 

insets), where m m = 0, are plotted against L~ 1 '" with u^ 1 = 0.7. 
The values of the cumulant crossing points in (a), (b), (c) are taken 
to correspond the respective transition points for L — > oo. 

It might be noted that in the low temperature and At 
regions, there seems to be significant finite size scaling of 
the transition (m w = 0) point (see the insets of Fig. 2a, 
b, c). In fact, in Fig. 2d, the finite-size scaling analysis of 
those data is presented. For the other cases, there seems 
to be no significant finite size effect on the transition 
point (cf. insets of Fig. 3a, b, c), indicative of a mean- 
field nature of the transition in this range. It may be 
noted that to compare the finite size effects, we normalise 
the parameters T or h p by their ranges required for full 
magnetisation reversal. In fact, this weak finite size effect 
for high T and At regions did not lead to any reasonable 
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FIG. 3. Behavior of U near the transition, driven by (a) h v at 
a fixed value of T (=2.0) and At (=5) with sequential updating, 
(b) T at a fixed value of h p (=0.5) and At (=10) with sequential 
updating, and (c) h p at a fixed value of T (=1.5) and At (=5) with 
random updating, for different L, averaged over 1000 to 6000 initial 
configurations. The fluctuations are smaller than the symbol size. 
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The insets show the typical behavior of the magnetisation m w at 
the time of withdrawal of the field pulse by varying (a) h p at a 
fixed T and At, for L = 50 and 400, (b) T at a fixed h p and At, 
for L = 50 and 200, (c) h p at a fixed T and At, for L = 50 and 
200; m m = at the transition point. 

For the static transition of the pure two-dimensional 
Ising system, U* ~ 0.6107 |7J,|8|,|!|- For low temperature 
(and low At) regions of the magnetisation-reversal phase 
boundary, the observed values of U* (in the range 0.42 
- 0.46) are considerably lower than the above mentioned 
value for the static transition. There is not enough in- 
dication of finite-size effect in the U* value either (cf. 
0)- This suggests a new universality class in this range. 
Also, the finite-size scaling study for the effective transi- 
tion points here (see Fig. 2d) gives a correlation length 
exponent value (y ~ 1.4) larger than that of the static 
transition. For comparatively higher temperatures (and 
high At), the U* ~ 0+ at the crossing point. Such small 
value of the cumulant at the crossing point can hardly be 
imagined to be a finite-size effect; it seems unlikely that 
one would get here also the same universality class and 
U* value will eventually shoot up to U* ~ 0.44 (for larger 
system sizes), as for the other range of the transition. On 
the other hand, such low value of U* might indicate a 
very weak singularity, as indicated by the mean field cal- 
culations mentioned in the introduction. In fact, even 
for the static transition, as the dimensionality increases, 
and the singularity becomes weaker (converging to mean 
field exponents) with increasing lattice dimension, the cu- 
mulant crossing point U* decreases (U* ~ 0.61 in d = 2 
to U* ~ 0.44 in d = 4) @. We believe the mean field 
transition behavior here, as mentioned earlier, is even 
weaker in this dynamic case as reflected by the value 
U* ~ 0+, corresponding to a logarithmic singularity (as 
in eqns. (J5J and 0). 

IV. SUMMARY AND CONCLUSIONS 

The universality class of the dynamic magnetisation- 
reversal transition, induced by a competing pulse, in an 
Ising model on a square lattice, below its static ordering 
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temperature, is studied here using Monte Carlo simula- 
tions. Both sequential and random updating have been 
used. The phase boundary at any T (< T°) is obtained 
first in the h p — At plane. They of course depend on 
the updating rule. The phase boundaries obtained com- 
pare well with the nucleation theory estimate /i p lnAt = 
constant along the boundary The mean-field the- 
ory applications 0, 0| indicated time and length (eqns. 
and respectively) scale divergences at these phase 
boundaries. Extensive Monte-Carlo studies for the fluc- 
tuations in the order parameter \m w \ and internal ener- 
gies etc. showed prominent divergences along the phase 
boundaries 0. Fourth order cumulant (U) of the order 
parameter distribution is studied here for different sys- 
tem sizes (upto L = 800) around the phase boundary 
region. The crossing point of the cumulant (for differ- 
ent system sizes) gives the transition point and the value 
U* of the cumulant at the transition point indicates the 
universality class of the transition. In the low temper- 
ature and low pulse width range, the U* value is found 
to be around 0.44 (see Figs. 2a, b, c). The prominent 
discripancy with the U* value (~ 0.61) for the static 
transition in the same model in two dimensions indicates 
a new universality class for this dynamic transition. In- 
deed, the finite-size scaling analysis (Fig. 2d) suggests a 
different (larger) value of the correlation length exponent 
also. For comparatively higher temperatures and higher 
pulse widths, the U* values are very close to zero (see 
Fig. 3a, b, c), and the transitions here seem to fall in 
a mean-field-like weak- singularity universality class sim- 
ilar to that obtained earlier and indicated by eqns. 
(0 and J7J. Here, the finite size effects in the order pa- 
rameter and the transition point are also observed to be 
comparatively weaker (see insets of Fig. 3). 
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